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' The passage of light or of electrons through a disordered medium is modified 

(N 

CN ' in the presence of resonances. We describe a simple model for this problem, 

and present first results. 
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During the last decade, many phenomena related to the passage of waves through a 
disordered medium have been intensely studied. This holds both for electrons, i.e. for the 

g 

o . 

waves. As examples, we mention conductance fluctuations in mesoscopic probes, and weak 



amplitude waves of the Schrodinger equation, and for light, i.e. for classical electromagnetic 



and strong localization effects including the enhanced backscattering of light reflected from a 
disordered medium. For electrons, disorder is typically caused by scattering from impurities. 
For light, disorder is often produced artificially: A powder is immersed into some liquid. 
Random scattering occurs if the indices of refraction of powder and liquid are sufficiently 
different. 

The present paper has been motivated by an experiment performed several years ago 
||T|] on the scattering of light by a disordered medium. In this experiment, the powder used 
(TiC^) consisted of grains with a size distribution centered around a diameter of 220 nm. 
For such granules, a Mie resonance occurs close to the wavelength A ~ 630 nm of the laser 
light used in the experiment. This led to a resonance enhancement of the scattering. The 
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effects of such enhancement were seen by comparing the diffusion constant D (determined 
from the intensity autocorrelation function versus frequency of the transmitted light), and 
the transport mean free path I (determined either from weak localization, i.e. from enhanced 
backscattering, or from the dependence of the transmission intensity on the length L of the 
disordered slab). One expects D and / to obey the relation D = 1/3 veI, with ve the energy 
transport velocity through the disordered medium. In the present case, the ratio 3D /I 
yielded a value v e = (5 ± 1) 10 7 m s _1 for the transport velocity ve which is about an order 
of magnitude smaller than the phase velocity. This surprising result has been understood 
both for small concentration |J and, more recently, also for strong concentration of the 
scatterers. Qualitatively speaking, the transport velocity is reduced because on its way 
through the medium, the energy is stored for some time in the Mie resonances. 

In this paper, we study the passage of waves through a medium with randomly positioned 
resonant scatterers from the point of view of random matrix theory (RMT). To the best of 
our knowledge, this has not been done before. Our motivation is the following. First, the 
approach taken in refs. JIJ and is not entirely systematic. Indeed, it combines results of 
diagrammatic impurity perturbation theory (used for analysing the intensity autocorrelation 
function) with the use of the Bethe-Salpeter equation |IJ , or of a mean-field approximation 
||, for the calculation of the transport velocity ve- It is our hope that a suitable modifica- 
tion of RMT (capable of accounting for the presence of resonance scattering) may yield a 
uniform account for both, average and fluctuation properties. This is all the more desirable 
since during the last decade, RMT has developed into one of the most important tools for 
investigating wave propagation in disordered media. The problem studied in this paper calls 
for an extension of existing RMT and would thus widen the scope of this theory. Second, 
RMT would make it conversely possible to put the present problem into a wider context, 
and to study, for instance, the effect of resonances on localization properties. 

We have embarked on this project only recently. The present paper is, therefore, no 
more than a work-in-progress report. We simplify the problem by focussing on scalar waves 
with a simple dispersion relation (i.e., electrons), thus postponing the description of the 
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propagation of light waves to a later stage. The RMT model with resonances is introduced 
in section [II]. In section JTJ, Efetov's supersymmetric generating functionals || in the version 
of ref. Q are used to handle the technical problems resulting from the model. First results 
obtained in the simplest case are presented in section [TV]. 



II. THE MODEL 

To motivate the model, and to identify the relevant parameters, we begin with a few 
qualitative considerations. We consider a volume V of liquid containing a volume concen- 
tration p (defined as volume of grains per total volume) of resonant powder. To simplify the 
problem, we assume that all grains (i.e., all Mie scatteres) are exactly equal, each having 
volume $, and resonance frequency u)\. Then, the number m of grains in V is m = pV/Q. 
The total volume m$ = pV occupied by the grains is denoted by V gra i n , the remaining 
volume V rem left for the liquid is V rem = V(l — p). 

We denote by (ou) the density of states per unit volume in the liquid for a classical 
scalar wave of frequency u. We have || 

where vq is the phase velocity of the wave in the liquid. For p = 0, the total density of states 
in the liquid is pi = p^\oj)V. For p ^ we only have the remaining volume V rem at our 
disposal, so that pL,rem = PiVrem — PlV(1 ~ p) ■ F° r later use, we note that for a vessel of 
length L — 17 pm and area A = 6 10 3 (p m) 2 , for vq = c (the velocity of light in vacuum), 
and for a wave length A = 630 nm, the mean level spacing = p~j} is about 4 10~ 7 eV. 

With increasing concentration of scatterers, the remaining volume shrinks, and with it 
the level density of the liquid. We may ask: Where is the "missing level density" p° L V pi To 
answer this question, we consider a single grain. In this grain, there are many Mie resonances, 
with a level density pMie which is proportional to <l> (so that pMie = P°Mie® ^ Pl)- The total 
level density in all m grains is mp Mie = p° Mie m& = p° Mie Vp. This is exactly the missing 



3 



contribution provided we have p° L = p% ie - In practice, this condition will not hold, because 
the level density depends on the index of refraction of the medium. For purposes of the 
present qualitative discussion, we disregard this fact. 

We see that the immersion of the grains in the liquid causes a concentration of level 
density at the position of the Mie resonances: While p L is lowered from p° L V to p° L V(l — p), 
PMie increases from zero to p° Mie Vp. We note that on the scale of the mean level spacing 
oIl = 1/ (p° L V(l — p)) in the liquid, puie is not smooth but sharply peaked at the position of 
the Mie resonances. The area under each peak is m. We also note that typically, m 3> 1: 
For a concentration of 1 percent, a grain size of 200 nm, and a volume of 10~ 13 m 3 , we have 
m ~ 10 5 ! 

We expect the shape of the peak of pMie to be Lorentzian, with a width Y given by the 
decay width of the individual Mie resonances (we assume that the concentration is sufficiently 
low so that the resonances do not "talk" to each other). Noticeable optical effects can be 
expected only if the peak height 2m/{iiT) is at least comparable with p L , i.e. for mdi ~ T, 
or for p > pi^T . When is this condition met? For light, we typically have V ~ 0.1 eV. Using 
this figure and the values for pi and $ given above, we find p > 2 percent. 

Naturally, these considerations are only semiquantitative. Nonetheless, they establish the 
following picture. For grain concentrations in the percent region, there occurs a significant 
local enhancement of the level density at the characteristic frequency uj\ of the Mie scatterers, 
with a width given by the decay width T of the Mie resonance. Typically, this width is very 
much larger than the mean level spacing di in the liquid. For frequencies u differing from 
the Mie resonance frequency lq\ by more than a few units of T/h, the wave incident on a 
grain does not excite the Mie resonance, and the scattering by the wave on the grains is a 
totally random process kin to diffusive impurity scattering. As uj approaches u>i, resonant 
scattering will be superposed onto this random process. Even for uj — u>i, there will still 
be random scattering since upon hitting a Mie scatterer, only part of the incident wave will 
excite the resonance, the remainder undergoing hard-sphere reflection. 

How can we model this situation in the framework of RMT? We recall how diffusive 
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scattering in the absence of Mie resonances is modelled in the framework of RMT, for a 
quasi one-dimensional system of length L ||: The system is considered as consisting of 
many longitudinal slices of length I <C L. Within each slice, the Hamiltonian is modelled as 
a random Hamiltonian matrix (a member of the random-matrix ensemble with the proper 
symmetry). Neighboring slices are coupled by Gaussian-distributed uncorrelated random 
matrix elements. The first and the last slice are coupled to the leads or, equivalently, to 
the channels. Elements of the scattering matrix connecting incident and outgoing channels 
define the conductance, i.e. the transmission through the disordered region. 

With proper modifications, we use this model also for the present problem. To begin 
with, we consider only Schrodinger waves. The modifications caused by the different form 
of the Hamiltonian for light waves will be addressed at a later stage. We retain the idea 
of a division of the quasi one-dimensional system into slices. To account for the presence 
of Mie scatterers, the form of the Hamiltonian has to be modified: Within each slice, our 
Hamiltonian H is a matrix of dimension N + ms where ms is the number of Mie scatterers 
within each slice. The Hamiltonian H contains a purely random part (a matrix Hqoe of 
dimension N belonging to the Gaussian orthogonal ensemble [GOE]), a diagonal part of 
dimension mj with diagonal elements E% = %uj\ corresponding to the presence of mj Mie 
scatterers with equal resonance frequencies u\, and a rectangular matrix V which couples 
the ms resonaces to Hqoe- Since we will focus on the Hamiltonian for a single slice in all 
that follows, we replace mg by m for simplicity. With I m the m-dimensional unit matrix 
and with T denoting the transpose, we thus have 



H 



( Hgoe V ^ 



y V T E 1 XI m j 



(2) 



The dimension N of Hgoe is taken to infinity at the end of the calculation. The vari- 
ance of the matrix elements Hgoe is normalized as usual: The ensemble average of 
(H G oe)^u{Hgoe)^'u' is given by (X 2 / N){5^5 u y + S^yd^y). Here, A (not to be confused 
with the wave length of the scattered light) has the dimension of an energy and determines 
the average level spacing ds in each slice. For m — 0, the average level density obeys the 



semicircle law, and the dependence of ds on energy E is given by n\/Nyj 1 — (E/(2X)) 2 . For 
m ^ 0, we expect the average level density to deviate from the semicircle law by aquiring 
an additional peak at energy E\. The width of the peak will be determined by the matrix 
V. By assuming the m-dimensional matrix in the lower right-hand block to be diagonal, we 
implement the assumption that the m resonances do not "talk" to each other, i.e., are not 
coupled directly. By the same token, we assume that resonances in different slices are not 
directly coupled to each other. Both assumptions are obviously valid only for a sufficiently 
low concentration of Mie scatterers. Because of this assumption, the rectangular matrix V 
can, without loss of generality, be taken to be diagonal. Indeed, V can be written in the form 
0\Vr>0\ where Vb is diagonal, and where 0\ and O2 are orthogonal matrices of dimension 
N and m, respectively. Transforming H with the orthogonal matrix O = diag(Oi, O2), and 
absorbing the matrix 0\ into Hgoe by virtue of the orthogonal invariance of the GOE, we 
obtain a new Hamiltonian of the form of Eq. (fj), but with V replaced by the diagonal matrix 
Vjj- Since we assume all scatterers to be identical, we take Vp to be a multiple of the unit 
matrix, Vd — v x 7 m . This completes the definition of the model. 



III. SUPERSYMMETRIC FORMALISM 

So far, we have only considered the simplest apects of the model defined in section ||: 
For technical reasons, we have replaced the GOE by the unitary ensemble (GUE), and we 
have calculated the mean level density and the two-point level correlation function. This 
was done to make sure that the model behaves as we expect, and to gain experience with 
a Hamiltonian of the form of Eq. (|2]). We have used the supersymmetry technique of refs. 

We omit all details, use the notation of ref. j|, and give only a few results which 
display characteristic differences between the usual procedure, and the present model. 

After averaging and the Hubbard-Stratonovich transformation, the exponent of the gen- 
erating functional for the two-level correlation function (with energies E ± u) has the form 
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E - Xa - loA - J -V 

E - Xa -uA- J 

-V E - E x - cuA - J 



}m 

}N - m • 

}m 

(3) 



The symbols J stand for the source terms. The dimensions of the block matrices are indicated 
at the side, and V stands for v x I m . Variation with respect to a leads to the saddle-point 
equation 

1 m E — Ei 



a 
A 



1 



m, 



+ 



(4) 



N'E-Xa N (E - Xa)(E - Ex) - v 2 ' 
As usual, we have assumed that uj ~ 1/iV, and we have omitted the (Infinitesimal) source 
terms. We take the limit N — > oo and keep m fixed. This yields the standard saddle-point 
equation 

1 



a 
A 



E — Xa 



(5) 



with solution 



E -a 
a =2X +iA > 



(6) 



where A = J 1 — (E/2X) 2 . In contrast to the usual case, the average level density is no 
longer given by the saddle-point solution of Eq. (|6]) but attains an extra piece which is due 
to the last block in the matrix displayed in Eq. (||). Likewise, the two-point function differs 
from the regular form, and has to be calculated by keeping carefully track of the extra terms 
appearing in Eq. (|3|). The results follows. 

For the average density of states, we define the level density for the Hamiltonian H as 
p(E) = J2i.d{E — Ei). Let a bar denote the ensemble average. Then, we find 

TA(E)/n 



-tt^ A(E)N 
p(E) = \ ' +m- 



(7) 



Att ' " (E - E x - (E/2X)T) 2 + (TA) 2 ' 
where T = v 2 / X defines the effective width TA of the Lorentzian. The form of Eq. (|7]) is 
expected. Indeed, 7rA/(A^A) is the mean level spacing d of the GUE and is to be identified 



with d,L, the mean level spacing in the liquid. The last term in Eq. (0) has area m, peak 
height m/(7iTA), and Lorentzian shape, all as anticipated in section [Q[ The occurrence of 
the factor A reflects the semicircle law for the GUE level density, and the factor E/(2X) 
multiplying T is an associated level shift which vanishes in the center of the semicircle where 
E = 0. The form of T is also in keeping with expectations. Indeed, we recall that the 
GUE matrix elements have variance X 2 /N. We expect the matrix elements v occurring in 
Eq. (U) to scale similarly with N, so that v 2 = a\ 2 /N where a is independent of N. Then, 
T = aX/N so that T is given as «7r/A times the mean level spacing of the GUE. For light 
scattering, a ^> 1 although we also consider the case a ~ 1 in this paper. 

The Lorentzian in Eq. (|7|) is plainly a superposition of m individual resonances. It may 
come as a surprise that the indirect interaction between resonances caused by their joint 
coupling to the random Hamiltonian does not cause any level repulsion between the m 
Lorentzians. To clarify this issue, we have solved the saddle-point equation (f|) to next 
order in iV _1 . We find that this leads to an energy-dependent broadening of the Lorentzian 
in Eq. (|7]) which obviously is of order iV -1 relative to the terms we have kept. We conclude 
that within the approximations leading to Eq. (|7|), the indirect coupling between resonances 
is negligible. 

Since the presence of the m resonances at energy E\ destroys translational symmetry, 
the evaluation of the level-level correlation function is slightly more complicated than is the 
case usually. We calculate 

C(E, u) = p{E + u)p(E -u)- p[E + u) p{E - u). (8) 

The disconnected terms cancel as usually if the "boundary terms" or "Efetov-Wegner" 
terms are taken into account. The remaining integrand contains an exponential with 
argument 

- 2i7ru/d trg(l + oti) + m trg log ((E - E l - ET/(2X)) 2 + A 2 T 2 - u 2 - 2^rA(l + e*i)) . 

(9) 
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Here, ot\ is defined as in ref. M. It is easily seen that for uj <C T, this expression reduces 



to —2i7ruj/d e ff(E) trg(l + cti) where d e ff(E) is the inverse of the average level density p(E) 
given in Eq. (|7|). By virtue of the integration over the variables in a±, the cu-dependence of 
C(E,u) is effectively cut off at u ~ d e ff. For values of V 3> d e ff, the normalized two-point 
autocorrelation function therefore essentially coincides with the one for the pure GUE except 
for a rescaling of the average level spacing. Differences can occur only in the wings of the 
distribution, or in cases where V ~ d e /f. The figures presented in the next section are based 
on the exact form of C (E, uS) given by 

1 / iD ^ - m+1 



C(E, u) = p(E + uj)p(E -uj) = Re 



exp(2irujD/d)T(-m + 1, 2ttuj(D - i)/d) 



2d 2 \2nu} , 

r 1 m / iD \ ~ m /D I \ 

x / d\ 2 (iD + A 2 ) m ex V (-2i7tu\ 2 /d) --—( — ) exp(2nuD/d) - + — 

J-i 2ird\2irujJ \uj 1 Ay 

/l rrP' f iD \ 

d\ 2 (iD + X^- 1 exp(-2t7TLoX 2 /d) - -— 
-1 67T Z \2tIUJJ 

xexp(2nujD/d)T(-m-l,2TTUj(D-i)/d) ^ d\ 2 {iD + A 2 ) m " 2 exp{-2muo\ 2 /d) 



i2 / „' n \ —m—l n2 



D 2 + 1 



^ 2 



(10) 



where T(a, b) is the Gamma function, and where 

D = [{E - E x - (E/2\)T) 2 + (rA) 2 - uo 2 }/{2uoT^). (11) 

The remaining integrals in Eq. (|T0|) can be calculated analytically for every m. The resulting 
expressions are too lengthy to be given here. 

In the sequel, we also consider the normalized quantity 

C*or{E,u>)- . (12) 



IV. RESULTS 

We have commented above on the form of the average level density given in Eq. ([7|); 
these comments will not be repeated here. Rather, we focus attention on the normalized 
level-level correlation function defined in Eq. (|T2"D. With all energies scaled in units of the 
GUE mean level spacing d = ttX/N at the center of the unit circle, we present in Figure 1 
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the function C nor versus u/d for parameter values as given in the caption. For comparison, 
we also show the standard GUE level-level correlation function sin 2 '(2ttuj / 'd e ff) / '(2-nuo / 'c? e //) 2 • 
We note that in the figure, we have chosen Y = 0.3d so that differences between the two 
functions should be visible in the wings. As anticipated above, we find that the two functions 
nearly coincide within the peak. While the GUE correlation function is positive semidefinite 
everywhere, the function C nor oscillates in the wings around zero. 

In Figure 2, we show C nor for the case T = 3d. Here, there is no dis- 
cernible difference between this function and the standard GUE autocorrelation function 
sin 2 (2ttu I d e f /) /(2ttuj/ d e ff) 2 provided that the latter is scaled with d e ff. The second func- 
tion in Figure 2 with the wide peak shows the GUE correlation function evaluated without 
rescaling, i.e. the function sin 2 (27ra;/(i/(27ra;/(i) 2 . The difference between both functions is 
substantial. 

In summary, we have presented first analytical results for a random-matrix model de- 
scribing wave propagation in a random medium with resonant scattering. We are confident 
that this work will produce further interesting results. 
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Fig.l. The normalized level-level correlation function of Eq. (8) for m = 5, T = 0.3d, 
E = E x = (full line) and sm 2 (2irLu/d eff )/(2iiLu/d eff ) 2 with d eff = d / {I + m / (i\T)) (dotted 
line). 
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Fig. 2. The normalized level-level correlation function of Eq. (8) for m = 5, T = 3d, 
E = Ei = (which coincides almost exactly with sm 2 (2nuj/d e ff)/(27iuj/d e ff) 2 with 
d e ff = d/(l + ra/(7iT))) (full line) and the standard normalized GUE correlation function 
sin {2-kuj / d) / {2nu / d) 2 (dotted line). 
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